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We study periodically driven bosonic scalar field theories in the infinite N limit. It is well-known 
that the free theory can undergo parametric resonance under monochromatic modulation of the 
mass term and thereby absorb energy indefinitely. Interactions in the infinite N limit terminate 
this increase for any choice of the UV cutoff and driving frequency. The steady state has non-trivial 
correlations and is synchronized with the drive. The 0{N) model at infinite N provides the first 
example of a clean interacting quantum system that does not heat to infinite temperature at any 
drive frequency. 

PACS numbers: 


I. INTRODUCTION 

Experiments in cold atomic systems have generated 
much interest in the dynamics of periodically driven 
many-body Hamiltonians mm- Energy is not conserved 
in such systems; instead by the Bloch-Floquet theorem, 
the eigenstates have the form: 

m)) = e-^^^\m ( 1 ) 

where e is a quasi-energy defined modulo the fundamen¬ 
tal drive frequency and \4>{t)) has the same periodicity as 
the drive [3]. There are two intimately related questions 
in such a system: i) is there a late time steady state with 
non-trivial correlations and finite energy density, and ii) 
is such a late-time steady state synchronized with the 
drive [37!? Already for a single two-level system, the an¬ 
swers to these two questions are non-trivial as the system 
can coherently Rabi flip-flop at a different frequency from 
that of the drive. 

At the next level of complexity are integrable many- 
body systems such as periodically driven non-interacting 
fermions. The dynamics is governed by an effective 
quadratic Floquet Hamiltonian; thus the stationary state 
coincides with an appropriate periodic generalized Gibbs 
ensemble There has been much theoretical progress 
classifying the topological structure of Floquet bands ni¬ 
ls] and experimental progress studying such states in cold 
atomic systems |5| and topological insulators nnuni. 

The nature of the steady state in driven interacting 
systems is less clear. Standard linear response theory 
suggests that any finite frequency drive heats the system 
to infinite temperature. That is, the local reduced den¬ 
sity matrix approaches the identity. Refs. HMD argue 
for this scenario in generic ergodic systems with locally 
bounded Hilbert spaces. These findings are in contradic¬ 
tion with Refs. nMii that claim that certain spin mod¬ 
els do not heat to infinite temperature when the drive 
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frequency is above a finite threshold. For strongly disor¬ 
dered spin systems whose time independent Hamiltonian 
is many-body localized, several recent studies find the 
same threshold behavior [TBl - ED] . For systems with lo¬ 
cally unbounded Hilbert spaces, even less is known. One 
recent study of the many-body Kapitza pendulum hnds 
threshold behavior |21j : we compare their results to ours 
later. 


In this article, we approach this problem using the 
large N expansion for interacting bosonic systems. The 
0{N) model at infinite A^ is a canonical model for 
symmetry-breaking in statistical mechanics [22j . Its equi¬ 
librium properties are exactly soluble and capture the 
correct topology of phase diagrams in various dimensions. 
It is also a canonical model for the unitary dynamics of 
interacting theories and a workhorse of many fields in¬ 
cluding cosmology and condensed matter [23)130] . The 
Floquet dynamics in the infinite N limit is the focus of 
this work. We comment on 1/A^ corrections towards the 
end, but reserve a full-blown treatment for the future. 


We begin in Sec. ]TT] by reviewing the response of 
the Gaussian model under periodic driving of the bare 
mass in the paramagnetic phase, previously studied in 
Refs. [3T| [02] . In this theory, the energy density grows 
exponentially in time via the parametric resonance of se¬ 
lected momentum modes [33j . Incorporating the effect of 
interactions at A^ = oo removes this divergence for any 
fixed cutoff and drive frequency (Sec. HI). The long-time 
steady state synchronizes with the drive, but depends on 
initial conditions because of the integrability of the the¬ 
ory. In Sec. EYi we turn to the symmetry-broken phase. 
Again, interactions prevent indefinite heating, although 
the time averaged magnetization decays to zero. Oddly, 
the late time magnetization can oscillate at half the fre¬ 
quency of the drive; we argue that this is an artifact of 
the infinite N limit. The 0{N) model at infinite N is the 
first example of a many-body system that fails to heat to 
infinite temperature for any periodic drive. 
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FIG. 1: Left: the energy spectrum of a Id particle in a cosine 
potential vs quasi-momentum q. The shaded region is the 
mode range. If the mode range includes band-gap(s), the 
Gaussian model is unstable. Right: the effective Mathieu 
spectrum in the driven 0{N) model at late times. The mode 
range always lies within a band. 


II. GAUSSIAN MODEL 


and n(fc,i) imposes the following constraint: 

Im[4^(t)/;(t)] = l/2 Wk (4) 

The mode functions satisfy the equations of motion: 


ddP' 


+ \k\ -b To - ri cos( 7 t) /c = 0 


(5) 


At each momentum k, this is the well-known equation of 
motion of a parametrically driven two-dimensional har¬ 
monic oscillator (as /g is complex). The transformation, 
t 2 t/ 7 , —>• \/ 2 / 7 /g, makes time and the mode func¬ 
tion dimensionless and brings Eq. ([^ to the canonical 
form: 


^+Ek-2g cos( 2 t) ) /j = 0 


El. = 


4(|fcp+ro) 


2ri 

9 = ^ 
T 


( 6 ) 

(7) 


This differential equation is the Mathieu equation |55] . 

The Mathieu equation is a familiar beast in band the¬ 
ory; it is the Schrddinger equation of a one-dimensional 
particle in a cosine potential. This identification provides 
a dictionary between the mode functions f^{t) and the 
Bloch wavefunctions of the Schrodinger equation: 


Consider the free 0{N) model in d dimensions with a 
sinusoidally varying bare mass. The Hamiltonian is: 


Ho{t) = 


rA ,d 


d‘^k 


|n.(fc)|^ 


+ (|fcp+r(t)) 




r{t) = tq — ri cos( 7 t) 


( 2 ) 


where ^i{k) and nj(fc) satisfy the canonical commutation 
relations: 


[$,(fc),n,(fcO] = i{2ny6{k-E)S,j. (3) 

The component index i runs from 1 to the ultraviolet 
cut is given by A in momentum space and the bare mass 
is r{t). We suppress the label i as the initial states are 
0(iV)-symmetric. In equilibrium in the absence of drive 
(ri = 0 ), the model is paramagnetic for rg > 0 at all 
temperatures and ill-defined for rg < 0 . 

In the bulk of the article, we analyze the response un¬ 
der a monochromatic drive at a single frequency 7 . As 
the Gaussian theory is a linear system, the response to a 
polychromatic drive follows by superposition. The inter¬ 
acting 0{N) model is however not a linear system. We 
address generic periodic drives in this model at the end 
of Sec. mil 

In the Gaussian model, the equations of motion for 
t), n(fc, t) are linear. We expand these operators in 
a fixed basis of creation and annihilation operators, e.g. 
^{k,t) = This defines the complex 

mode functions The commutation relation between 


t —)■ Spatial coordinate 
—)■ Wavefunction 

Ek —>■ Energy of the wavefunction ( 8 ) 

A number of properties of the Mathieu spectrum (Fig.[^ 
follow directly from this mapping: 

1. By the Bloch theorem, the spectrum of the par¬ 
ticle is labelled by a quasi-momentum q G (— 1 , 1 ] 
and a band index m = 1 , 2 ... with eigenfunctions: 

= e“*'^“(/)'"(a:) where = (j/^[x + n). 

2. When the amplitude of the drive g is zero, the spec¬ 
trum is degenerate at the centre and the edges of 
the Brillioun zone. Any g ^ 0 opens a gap at these 
degenerate points; the mth gap is approximately 
g'^/{{m/2)\)^ for large m. The gaps decrease faster 
than exponentially in m. 

3. The bandwidth of the mth band is approximately 
2m -b 1 at large m. 

Using the mapping in Eq. the fc-modes in the 
Gaussian theory sample energies from Eq = drg/y^ to 
E\ = 4(rg -b A^)/ 7 ^ in the Mathieu spectrum (shaded 
region in Fig. [^. We call the range of energies between 
Eq and Ea as the ‘mode range’. If the mode range lies 
within a Mathieu band, then by (1) above, each mode 
function is a superposition of the two solutions at Eqk 
and is oscillatory in time: 

Vr (0 + (9) 
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where ak , /3fc are complex numbers determined by the 
initial conditions and Eq. Q. The energy density and 
other spatially local observables involve integrals over the 
mode functions in fc-space. As the magnitude of each 
mode function is bounded in Eq. ([^, all such observables 
remain bounded as t —>■ (X). Further, it is straightforward 
to show that local observables oscillate in synchrony with 
the drive as t —)■ oo. For example, (d>f (<)) is given by: 

(‘J>-W) =f^ + 2Re[e-^*^'=‘xFW] 

where Sk and X^(^) are related to the parameters in 

Eq. The second term decays as as t ^ oo 

while the first term has the same period T = tt as the 
drive. Thus, ($|(t)) = + tt)) as t —>■ oo. 

The phase diagram is sketched in Fig. The simplest 
stable phase lies at small A^/ 7 ^,ro/ 7 ^; here the drive 
frequency 7 is much greater than any other energy scale 
in the system. However, the additional stable phases 
are non-trivial consequences of the band structure of the 
Mathieu equation. 

When the mode range intersects the Mathieu band- 
gaps, the corresponding mode functions increase expo¬ 
nentially in time, exhibiting parametric resonance. Con¬ 
sequently, the energy density (and other local observ¬ 
ables) also grow exponentially in time and the system 
heats indefinitely. The heating time-scale is given by 
the inverse of the largest band-gap intersecting with the 
mode range. It is important to note that at any fixed bare 
mass and drive parameters, the theory is always unstable 
for sufficiently large cutoff A as the mode range increases 
with A. 

It is sometimes useful to think of the evolution over 
a period T = 27r/7 as being generated by an effective 
Floquet Hamiltonian Hp\ 

U{T) = (10) 

where U{T) is the evolution operator for a period [51IM]. 
As the theory is quadratic, Hp can be chosen to be 
quadratic in the field operators. In the stable regime, 
the mode spectrum of Hp is non-negative and the eigen- 
modes are normalizable. By expanding any initial state 
in this eigenbasis, it is easily seen that the late time re¬ 
sponse is stable and periodic. In the unstable regime, 
on the other hand, the mode spectrum includes negative 
energies and corresponding unnormalizable eigenmodes. 
This is what allows the system to absorb energy indefi¬ 
nitely. For more details, see Ref. [32] . 

We end with three comments. First, the physics dis¬ 
cussed above applies to any spatial dimension d > 1. 
Second, in the phase diagram in Fig.[^ the stable region 
persists to some < 0. Thus, the driven Gaussian 

theory can be stable even when the equilibrium theory 
is not. Finally, “energy density” when unqualified refers 
to either the instantaneous energy density or the energy 
density with respect to the time-averaged Hamiltonian. 
Both diverge when the system heats up to infinite tem¬ 
perature. 
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FIG. 2: The dynamical phase diagram of the driven Gaus¬ 
sian theory at fixed g/ 7 ^. Unshaded/shaded regions are sta¬ 
ble/unstable. The band gaps in the Mathieu spectrum deter¬ 
mine the phase boundaries. 


III. DRIVEN PARAMAGNET 


We now turn to the interacting 0{N) model with a 
sinusoidally varying bare mass. The Hamiltonian reads: 

A ^ r 

H{t) = Ho{t) + — ^ / (<i>.(f))4 (11) 

At infinite N, ($i(T)$i(af))/V acts like a classical time- 
dependent field and can be replaced by its expectation 
value. In the absence of the drive in equilibrium (ri = 0), 
the model is paramagnetic for all rg > Tc and sponta¬ 
neously breaks the 0{N) symmetry for rg < rc- The 
value of Tc is determined by d: in d = 1 , Tc = — 00 , 
while in d > 2, rc is negative and finite. Further, the 
symmetry-broken phase extends to finite temperatures 
for d > 3. In this section, we focus on the coherently 
driven paramagnet. 

Expanding in a fixed basis of creation/annihilation op¬ 
erators as before, we obtain the equations of motion: 

( 12 ) 


assuming that fi^(t) is the same for every component and 
using the relation (no summation on i): 


mt)) 


r d<^k 

J (27r)^ 


m)\" 


(13) 


As the system is spatially homogenous, (<i)f) only de¬ 
pends on t. For more details, see Ref. [501 • As a conse¬ 
quence of the quartic term, Eq. (12) is non-linear in the 
mode-functions. Re-writing in dimensionless units: 



+ Ek 


2gcos{2t) + u 




hit) = 0 

(14) 
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where u = 8 A/ 7 ^ and Ek,g are defined in Eq. Q. 

Observe that m = 0 corresponds to the driven Gaussian 
model. When the Gaussian model exhibits stable behav¬ 
ior, it is clear that a small u merely dresses the steady 
state. The main result is that the stability persists at all 
parameters (even when the Gaussian model is unstable). 
The energy density always plateaus to a finite value at 
late times and the wavefunction has non-trivial correla¬ 
tions that can be described within an effective Gaussian 
model. Below, we present the intuition underlying this 
stability and numerical results that support this claim. 
We then construct explicit Floquet solutions to Eq. (14) 


at low drive frequency within the WKB approximation. 

The quartic term acts as a self-consistent correction 
to the energy E^. in Eq. (14). Define the instantaneous 
energy: 


Ek{t) — Ek + u($^(t)). 


(15) 


This identifies an instantaneous mode range. Suppose 
that at t = 0 , the mode range includes band-gaps in the 
Mathieu spectrum (Fig. [^. Then, the associated mode 
functions grow exponentially in time, and by Eq. (131, so 
does (<i)f(t)). Finally, this implies the mode range itself 
drifts up with time. As the bandwidth of each Mathieu 
band is proportional to its index, at large enough Eoit), 
the mode range lies within a single band (Fig. [^. We 
expect that all the mode functions then become oscilla¬ 
tory, the time averaged value of (4>^(t)) plateaus, and the 
system stops heating. Further, just as in the Gaussian 
model, we expect that all observables synchronize with 
the drive. 

A more refined version of the above argument corrects 
the drive at late times. At late times, (4>f(<)) is in syn¬ 
chrony with the drive: 




uj—0 




(16) 


where F[($|)](a;) denotes the Fourier amplitude at fre¬ 
quency ui. The mode functions thus satisfy a generalized 
Mathieu equation with parameters: 

Ek=Ek + uE[{<^>^M0), g = g - uE[{<^>^,)]{2) (17) 

The higher harmonics in ($f(t)) can be ignored as their 
amplitudes are small as compared to g. We expect that 
the mode range lies within a single band of the Mathieu 
spectrum with the parameters in Eq. 

The arguments above, while appealing, are not deci¬ 
sive. At short times, they assume that Eq. (14) can 


be treated as an effective Mathieu equation even though 
increases exponentially in time. At long times, 
they ignore the higher harmonics in the effective drive. 
We therefore turn to numerical simulations to confirm 
this picture. 

a. Numerics — Consider the d = 1 0{N) model in 
its paramagnetic ground state at t = 0. At time t = 0, 
we turn on the drive. With the parameters chosen in 


Fig.i the mode range at t = 0 includes the first band- 
gap. The modes near the center of the band gap set 
the initial time-scale of exponential growth of (<l>^(t)) to 
be approximately 2/g. The width of the mode range is 
1, while the bandwidth of the second Mathieu band is 
approximately 3; thus we expect (4>f (t)) to stop growing 
when Eq is approximately the energy at the bottom of the 
second band of the effective spectrum for the parameters 
in Eq. ( [T^ 

In FigT^ we plot (4)^(t)) vs the dimensionless time t 
at early (top) and late (middle) times. With g = 0.3, the 
initial time-scale of exponential growth is approximately 
7, in good agreement with the top panel in Fig. The 
growth rate is seen to decrease with time, as predicted 
by the effective Mathieu picture. At t/n « 10, (4>f(t)) 
saturates. This is evidence for stability; indeed, (4)|(t)) 
oscillates about a mean value for the next thousand peri¬ 
ods (middle panel). To connect to the effective Mathieu 
spectrum, we plot the Fourier spectrum (bottom panel). 
We see that i) (4’^(t)) is synchronized with the drive 
(Eq. @), ii) the harmonics at frequencies above two 
are suppressed, and iii) the value of Eq extracted from 
the figure coincides with the value of the lowest energy in 
the second band of the effective Mathieu spectrum. The 
energy density (not shown) is also finite and synchronized 
with the drive as t —>■ oo. Thus, the numerical simula¬ 
tions confirm the effective Mathieu picture and provide 
evidence for a steady state stabilized by interactions in 
the driven 0{N) model. 

b. Approximate Floquet solutions — We construct 


explicit normalizable solutions of Eq. (14) as the drive 


frequency 7 approaches zero within the WKB approxi¬ 
mation [35]. Physically, the system is most susceptible 
to indefinite heating at low drive frequency; the construc¬ 
tion of normalizable Floquet states in this limit is strong 
evidence that the model exhibits stable behavior at any 
drive parameters. 

Following the structure in Ref. |35j . we re-arrange 
Eq. ([I4|) as: 


Y d? 
4 dt"^ 


+ 7 (|A:p -f To - ri cos( 2 t)) -h 2X{^1[t)) ) f^{t) = 0 

(18) 


Note that we are still working in dimensionless units. To 
leading order in 7 in the WKB series, the mode functions 
are independent of k and are given by: 




1 


2V8A($f(t)) 


: exp 


±.1 


8A($f(t'))dt' 


where ($?(t)) is determined by the self-consistency con¬ 
dition in Eq. p|: 


{^m-K 


^2d/3 

Ai/3 


(19) 


AT is a dimension dependent constant. To next order in 7 , 
the mode functions depend on k and ($^(t)) is corrected 
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FIG. 3: Top and middle panels: vs t/iv at early and 

late times when the paramagnetic ground state is driven start¬ 
ing at t = 0. Bottom panel: Fourier spectrum of over 

a thousand periods of the drive. Parameters: d = 1, ro = 0.1, 
u = 1, A — 1, g = 0.3 and system size L = 100. The energy 
density shows similar behavior. 


by a 7 COs( 2 t) term. This expansion is systematic; we 
may go to as high an order in 7 as we desire and construct 
asymptotically accurate oscillatory solutions to Eq. (14|. 

c. Comments — We end with two comments. First, 
the effective Mathieu picture can be generalized to the 
case of a polychromatic periodic drive with fundamental 
frequency 7 . This is because the higher harmonics at 
frequencies imq, to > 2 do not change the qualitative 
picture of bands in Fig. As in the monochromatic 
case, the mode range drifts up in time until it lies within 
a single band of the effective Mathieu equation. Thus, 
the driven 0{N) model reaches an interaction stabilized 
synchronized steady state for any periodic drive. 

Second, the late time steady state is not described by 
a single temperature. Instead, it is described by an emer¬ 
gent periodic generalized Gibbs ensemble (PGGE) asso¬ 
ciated with the effective Gaussian model [3D]. The con¬ 
served quantities that feature in the PGGE are the mode 
occupations of the effective Gaussian model with param¬ 
eters in Eq. (17). Let denote the mode occupation 
of the oscillator at momentum k. At late times, 
commutes with the unitary evolution operator U(t) and 


is periodic with the same period as the drive. As these 
quantities are conserved only at late times, the PGGE is 
emergent. 


IV. DRIVEN FERROMAGNET 


Like the driven paramagnet, the driven ferromagnet 
does not heat to infinite temperature at any drive param¬ 
eters. Rather, it reaches a stable paramagnetic steady 
state in any d > 2. Unlike the driven paramagnet how¬ 
ever, the time period of the observables at late-time can 
be double that of the drive. We believe that this lack of 
synchrony is a consequence of the infinite N limit, as we 
discuss below. 

Let us be more precise. For tq < rc in d > 2, the 
0{N) symmetry is spontaneously broken in equilibrium. 
Let the symmetry be broken along the 1 direction in 
order parameter space so that ($i(a;)) is the non-zero 
uniform magnetization. The excitations along the 1 di¬ 
rection are massive. The remaining (N — 1) directions 
in order parameter space are soft and support Goldstone 
modes. The system is in the ground state at t = 0 when 
the mass drive is switched on. 

At infinite N, the equations of motion involve two clas¬ 
sical fields: ($i(t)) and The new classi¬ 

cal field defines the magnetization M{t) = (<i>i(t))/-\/]V. 
Defining mode functions for components 2,... N as be¬ 
fore and going to dimensionless units (t —)■ 21 / 7 , fj: —>■ 
y/2/jfj:, M{t) —)• y/2/jM{t)), the equations of motion 
are: 


“ 25C0s(2t) = 0 (20) 

+Eo- 2gcos{2t) +rf{t^M{t) = 0 ( 21 ) 

rf{t)=uj ( 22 ) 


where Ek,g are defined in Eq. ([^ and u = SA/q^ as in 
the previous section. The feedback term to the bare mass 
is denoted by rf{t) and involves an extra classical field 
as compared to Eq. (14). 

The intuitive argument for stability at late times is 
analogous to the one in Sec. |III| if we define the instanta¬ 
neous energy as Ek{t) = Ek + rf{t). At t = 0, the mode 
range extends from Eq{0) = 0 to Ea(0). Eq{0) is exactly 
zero as the Goldstone modes are massless. If the mode 
range lies within the lowest band, then the solution is 
stable as the mode functions are oscillatory (Sec. 0 . If 
however the mode range includes band-gaps, then rf{t) 
increases in time until the mode range fits in a single 
band of an effective Mathieu spectrum with parameters: 


Ek = Ek + F[rf]{0) g = g-F[rf]i2) (23) 

Again, F[rf]{uj) denotes the Fourier amplitude at fre¬ 
quency CJ. 
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The first consequence of the above picture is that Eq 
exactly coincides with the bottom of an effective Math- 
ieu band in the driven ferromagnet. In contrast, in the 
driven paramagnet, Eq could lie anywhere within an ef¬ 
fective band as long as the mode range fits in the band. 
A simple example is an initial state with Ao(0) = 2 and 
^'a(O) = 3 for the parameters chosen in Fig.[^ As the ex¬ 
ample suggests, the difference stems from the tunability 
of £'o(0) in the paramagnetic case. 

A second consequence is that the time period of late¬ 
time observables can be twice the period of the drive. 
This is because M (t) satisfies an effective Mathieu equa¬ 
tion with parameters Eq and g at late times and Eq coin¬ 
cides with the bottom of an effective band. If the bottom 
of the band corresponds to quasi-momentum q = 0 (mod 
2) in Fig.[^ then observables have the same period as the 
drive at late times. If instead the quasi-momentum is I 
(mod 2), the period is twice that of the drive. As the os¬ 
cillations themselves are likely an artifact of the infinite 
N approach, we do not expect this behavior at any finite 
N (see Sec. 0 . 

Finally, the time average of M{t) is zero at late times, 
as it is the solution of a Mathieu equation. Thus the 
steady state is paramagnetic. Interestingly, this is true 
in any dimension d > 2, even when there is a finite tem¬ 
perature ordered phase in equilibrium. 

d. Numerics — Numerical simulations confirm the 
effective picture discussed above. Consider the d = 2 
0{N) model in its ferromagnetic ground state at t = 0. 
At time t = 0, we turn on the drive. The parameters are 
chosen such that the initial mode range includes the first 
band gap in the Mathieu spectrum. 

Figs. Hand show M{t) and rf(t) in the time and fre¬ 
quency domain. As expected, M{t) decays to zero and 
rf{t) grows until Eq coincides with the bottom of the 
second band. As the bottom of the second band corre¬ 
sponds to quasi-momentum q = ±1, the period of M{t) 
is twice the period of the drive (Fig. [^bottom panel). 

The Fourier spectrum of rj (t) is also shown in the bot¬ 
tom panel of Fig. The dominant frequencies in the 
Fourier spectrum of rf(t) are seen to be at a; = 0 and 
oj = ±3. The small weight at w = ±3 opens mini-gaps in 
the effective Mathieu spectrum at q = ±1/2, as expected 
for a Gaussian model driven at frequency one [35] • 


V. DISCUSSION 

Intuitively, a Floquet system stops absorbing energy 
from a monochromatic drive when a fraction of its modes 
saturate, as in hole burning. This holds for fermions at 
the Gaussian level, as the fermionic modes do not in¬ 
teract, irrespective of whether the modes are spatially 
delocalized or localized. For bosons however, individual 
modes can absorb energy indefinitely by parametric res¬ 
onance at the Gaussian level. The reader might expect 
that going beyond the Gaussian level in either case leads 
to indefinite heating, as interactions allow for the ex- 
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FIG. 4: Top and middle panels: M{t) vs number of periods 
tj-n a.t early and late times when the ordered ground state is 
driven starting at t = 0. Bottom panel: Fourier spectrum of 
the time series between tj-n — 1200 and t/n = 1600. Parame¬ 
ters: d = 2, To = —0.4, u = 1, A = 1, g — 0.1 and system size 
L = 100. 


change of energy between modes. One way to cut off this 
heating is with sufficient quenched disorder, so that the 
modes that are nearby in energy are far away in space 
and are unable to exchange energy and thermalize. In 
this article, we discussed a different mechanism to cut off 
the heating, while still preserving enough symmetries to 
prevent full ergodicity. The clean driven 0{N) model at 
infinite N always has a finite energy density as t —?> oo, 
irrespective of the strength and fundamental frequency 
of the drive. 

At finite N, on the other hand, the 0{N) model is 
believed to thermalize. Previous work |36| argues that 
corrections up to order l/N"^ are required to see true 
thermalization in sudden quenches. This suggests that 
the driven 0{N) model could indefinitely heat once cor¬ 
rections to this order are included. However, the time 
scale for this heating would be parametrically large in 
N, so that the steady state discussed in this article would 
be observable up to this time. A detailed study of this 
question will be presented in a future work. 

Recently, Citro et. al |2T] studied the thermalization 
of another scalar field theory: the driven “many-body 
Kapitza pendulum” or the driven sine-Gordon model in 
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FIG. 5: Top and middle panels: r/(l) vs number of periods 
t/n at early and late times for the same parameters as Fig. 
is driven starting at t = 0. Bottom panel: Fourier spectrum 
of the time series between t/rr = 1200 and t/rr = 1600. The 
energy density shows a similar behavior. 


d = 1. They find that for parameters corresponding 
to both equilibrium phases—gapped and gapless—of the 
model, there is a critical frequency of the drive, below 
which the system heats indefinitely, but above which the 
heating stops. Their conclusions stem from various ap¬ 
proximate methods, of which one leads to equations that 
resemble the infinite N equations studied in this paper. 
Currently we do not understand the difference between 
their results in the gapped/paramagnetic phase and our 
results at low frequencies. We note that another of their 
methods—the application of the perturbative renormal¬ 
ization group to the Floquet Hamiltonian—is problem¬ 
atic as the stability of the Floquet Hamiltonian depends 
on the sign of the irrelevant terms. Understanding the 
precise connections between their results and ours will 
significantly clarify the steady state behavior of generic 
driven interacting bosonic systems. 
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